Effective temperature and compactivity of a lattice-gas under gravity 
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A distinctive feature of mean-field glassy dynamics is 
a peculiar violation of fluctuation-dissipation relations 
which leads to the definition of a time-scale dependent 
"effective temperature" and the possibility of con- 

structing a non-equilibrium thermodynamics of glasses 
and dense granular media Effective temperature 

also appears in athermal systems, where the high packing 
density regime is attained by compression or by using a 
confining potential. A particularly interesting situation, 
which is relevant to the study of granular materials, oc- 
curs when the confining force is gravity. In this case a non 
stationary inhomogeneous density profile generally arises 
and the notion of effective temperature may be then not 
well defined unless suitable conditions are verified. In 
this note we explore the possibility of defining a global 
effective temperature in an abstract model of granular 
material under gravity and weak thermal vibration, and 
its relation with the Edwards compactivity. 

The model consists of a gas of N particles on a body 
centred cubic (bcc) lattice where there can be at most 
one particle per site. There is no cohesion energy among 
particles and the Hamiltonian is 



N 



Ho^mg^ hi 



(1) 



where g is the gravity constant, hi is the height of the 
particle i, and m its mass. At each time step a parti- 
cle can move with probability p to a neighboring empty 
site if the particle has less than ly nearest neighbors be- 
fore and after it has moved ||]. Here p — min[l, x^^'*] 
where Ah = ±1 is the vertical displacement in the at- 
tempted elementary move and x — exp{—mg/ksT). 
We set mg/ks = 1 and ly = 4 throughout. At high 
enough packing density, dynamical models of this kind 
possess an extensive entropy of blocked states (defined as 
configurations in which any particle is unable to move) 
whose derivative is the so-called Edwards compactivity. 
For this reason such models exhibit a slow compaction 
dynamics reminiscent of dense granular matter |7||^. It 
was found in particular that during compaction a gener- 
alized fluctuation-dissipation relationship is obeyed 
giving a first evidence of an effective temperature in this 
regime. In Ref. the drift contribution to the longitudi- 
nal mean-square displacement was ignored p^ , leading 



to claim that "all measures of vertical correlation and re- 
sponse lead to the impossibility of defining effective tem- 
perature" and that "the vertical drift due to compaction 
leads to contradictory results" |11|. Here we show that 
there are no such contradictory results: in the slow com- 
paction regime the drift brings no qualitative change in 
the generalized fluctuation-dissipation relation found in 
Ref. 1^, and its effect is substantially negligible at high 
packing density. 

The fluctuation-dissipation properties can be charac- 
terized by applying a random perturbation to the system 
at times t> t^: 



N 



(2) 



where = ±1 independently for each particle, e is small 
enough to probe the linear response regime, and Q is the 
step function. The integrated response function is then 
defined as: 
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{t,t^)^-Y.{uAh,{t)) , 



(3) 



where Ahi{t) is the height difference between the per- 
turbed and unperturbed i particle at time t. The angular 
brackets denote the average over the thermal noise while 
the overline denotes the average over the random force. 
The 'mean-square displacement' between two configura- 
tions at time tw and i > tw is: 



1 ^ / 



h,{t^) + h{t^)-h{t) ), (4) 



where the drift motion is taken into account by the av- 
erage height: 
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Fig. |l| reports a parametric plot of mean-square displace- 
ment vs. response function at waiting time = 2^'^: 
it clearly shows that the presence of a slow longitudi- 
nal drift does not prevent the existence of generalized 
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fluctuation-dissipation relation. In particular, no qual- 
itative change occurs in the characteristic broken-line 
pattern when the drift term h(t„) — h{t) in Eq. (Q) is 
neglected, while appreciable quantitative deviations be- 
tween the two sets of data (with and without the drift 
term) only appear when the measurement time is quite 
long. Similar results were also found in the so-called 
FILG model under gravity [O. 
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FIG. 1. Parametric plot of mean-square displacement 
B{t,t^) vs. response function 2rx(t,tw)/e, during com- 
paction dynamics (circle symbols). The system is pre- 
pared in a random loose packed state with average density 
Prip ~ 0.707, and evolves under gravity and thermal vibra- 
tion with X = exp(— 1/r) = 0.2. The perturbation is turned 
on at the waiting time tw ~ "i^^ and measurements are carried 
out for times t in the range [tw, iw -I- 10^]. The slope of the 
dashed line is 0.23 (to be compared with 0.20 obtained by ne- 
glecting the drift term (x symbols)). The solid line with slope 
one is the equilibrium fluctuation-dissipation theorem, which 
is recovered by removing kinetic constraints (star symbols). 
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FIG. 2. Non-equilibrium fluctuation-dissipation relation 
in a compaction experiment as in FIG. 1, at different waiting 
times t„. The slope of the straight lines is 0.18, 0.20, 0.23, 
and 0.27, for increasing waiting time. 



The second result of the numerical compaction exper- 
iment is reported in Fig. it shows that the general- 
ized fluctuation-dissipation relation is obeyed at different 
waiting times. From the parametric plot of x(i,tw) vs. 
B{t,tv,) one can define a time scale dependent effective 
temperature by means of the relation: 



^dyn(^l^w) ^ 



2x(i,iw) ' 



(6) 



provided Tdyn is constant on that time scale, and where 
now it is understood a possible dependence of T^y^ on the 
density profile. Indeed, during compaction the system 
develops inhomogeneous density profiles |@,||, as it also 
happens after a sudden compression in zero gravity [p^ , 
and one may wonder about their influence on the effective 
temperature. In Ref. ||] the stationary density profile was 
interpreted as formed by two parts: a lower flat part at 
critical density pc — 0.84, and an upper equilibrium part 
in which kinetic constraints play no role. Fig. ^ shows the 
temporal evolution of the density profile corresponding to 
the above compaction experiment: one observes that the 
upper part of the bulk density profile increases faster than 
the lower one (even when the former become denser than 
the latter), and that the bulk profile is far from being 
flat (see inset of Fig. At late time the contribution 
of the top free interface is small for weak vibration, and 
- if sizeable - it would make higher the slope T/T^y^, of 
the fluctuation-dissipation plot (i.e. smaller the effective 
temperature) . While the contribution of particles at the 
bottom is negligible as they do not evolve at all. 
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FIG. 3. Temporal evolution of the density profile during 
compaction dynamics {x = 0.2, and time t = 2^°+* for k — 1 
to 4). Inset: bulk density proflle. 

That the generalized fluctuation-dissipation relation 
is not affected qualitatively by an inhomogeneous den- 
sity profile can be understood in terms of the mean- 
field dynamical model introduced in |l3| , and further 
generalized to nonzero gravity in In both cases 

the long-time relaxation of the local density factorizes: 
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Pc — p{z,t) = f{z)g{t). Since the mean square displace- 
ment can be written as B{z,t,t^) = F{z) G(t^t^) on 
long enough time interval t — t^, the violation factor en- 
tails two independent contributions: a purely geometric 
factor and a purely dynamic one. The latter contribu- 
tion is only responsible for the intrinsic violation of the 
fluctuation-dissipation relation. The notion of effective 
temperature therefore seems to be still reliable provided 
a geometric factor is taken into account (in this specific 
case the global geometric factor entering Eq. ^ would be 

= J /(z)dz/ / F{z)Az). Notice that for a purely flat 
density profile there is no difference between 'horizontal' 
and 'vertical' observables, (and = 1). 

The question that naturally arises is whether the longi- 
tudinal effective temperature can be interpreted in terms 
of the Edwards measure, which should now be obtained 
hy fixing the density profile of the experimental situation 
one wishes to reproduce [ p^ . The numerical implemen- 
tation of this strategy is however not straightforward. A 
more pragmatic approach consists in fixing a few feature 
of the profile (as suggested in this case by the expres- 
sion of f), such as average density, average slope and 
so on. This is quite similar in spirit to the construction 
of restricted Edwards measure which has been recently 
exploited in Ref . [^,^ , and generally improves the com- 
parison with numerical experiments. 



there is no dependence upon gravity for the bcc lattice. 
This is not however the most general case: interestingly, 
we found that for the simple cubic lattice, the anisotropy 
due to gravity, (say along the direction 001), suppresses 
the first order character of the phase transition present 
in the Edwards measure at density below « 0.7, |17|. We 
have then estimated T-Ed^ip) from the relation: 
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dp 



dp 



(7) 



where s{p) = —plogp—{l—p) log(l— p) is the equilibrium 
entropy |Q. The two effective temperatures, Tgdw and 
Tdyn, are shown in Fig. ^ at several densities correspond- 
ing to the average bulk density profiles of the compaction 
experiment. It is clear that a ponderable comparison 
is possible only when other features of the density pro- 
file (e.g. the average slope) are fixed, but much work is 
needed to test this point. 

In conclusion, we confirm the existence of a general- 
ized fluctuation-dissipation relation in an abstract model 
of dense granular matter which exhibits non-stationary 
inhomogeneous density profiles. The occurrence of a lon- 
gitudinal effective temperature in the slow compaction 
regime has been justified by a mean-field dynamical 
model and its relation with the Edwards compactivity 
has been discussed. 
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FIG. 4. Comparison between effective temperatures: Tdy„ 
is measured during the compaction experiment from the fluc- 
tuation-dissipation relation at time tw. lEdw is approximated 
through the Edwards measure with a homogeneous density p 
such that the average bulk density profile at time t„ of the 
compaction experiment is pav(iw) = P- 

As a preliminary attempt to relate the effective tem- 
perature Tdyn, to the inverse compactivity Tgdw, we have 
computed the Edwards entropy Sgdw (p) to the lowest ap- 
proximation, i.e. just by fixing a homogeneous density. 
In passing, the definition of blocked configuration in pres- 
ence of gravity requires here some care: if one assumes 
a flat profile, it is not clear why blocked configurations 
should depend upon gravity, like in the definition adopted 
in Ref. Q. We have therefore explicitly checked that 



J.J. Arenzon, A. Barrat and S. Franz are acknowledged 
for discussions and comments. 



[1] L.F. Cugliandolo, J. Kurchan and L. Peliti Phys. Rev. E 

55 3898 (1997). 
[2] S. Franz and M.A. Virasoro, J. Phys. A 33 891 (2000). 
[3] A. Barrat, J. Kurchan, V. Loreto and M. Sellitto, Phys. 

Rev. Lett. 85 5034 (2000). 
[4] S.F. Edwards, in Granular Matter: An Interdisciplinary 

Approach, A. Mehta ed. (Springer- Verlag, Berlin) 1994. 
[5] W. Kob and H.C. Andersen, Phys. Rev. E 48 4364 

(1993). 

[6] Notice that there is no purely horizontal motion in a bcc 
lattice: any elementary (one lattice spacing) displace- 
ment is either upward or downward. 

[7] M. Sellitto and J.J. Arenzon, Phys. Rev. E 62 7793 
(2000). 

[8] Y. Levin et al., Europhys. Lett. 55 767 (2001); J.J. Aren- 
zon et al., in preparation. 
[9] M. Sellitto, Phys. Rev. E 63 060301 (2001). 

[10] J.J. Arenzon, private communication (July, 2001). 

[11] A. Barrat, V. Colizza, and V. Loreto, Phys. Rev. E 66 
011310 (2002). 

[12] M.H. Vainstein and J.J. Arenzon, in preparation. 

[13] L. Peliti and M. Sellitto, J. Phys. IV (France) 8 Pr6 49 
(1998). The evolution of the density profile is detailed in: 
M. Sellitto, J. Phys. Condens. Matter 1455 (2002). 



3 



[14] This is indeed suggested by the functional nature of the 
nonzero gravity lattice gas entropy. For example, in the 
non-interacting case, Hamiltonian Eq. (1), the equilib- 
rium entropy is 



s 



W = - 2^ {p(^) log p{z) + [1 - p{z)] log [1 - p{z)]} 



From the minimization of the free energy one then gets 
mgz Ss 

[15] J. Berg, S. Franz and M. Sellitto, Euro. Phys. J. B 26 
349 (2002). 



[16] A. Lefevre, cond-mat/0202376 



[17] The occurrence of a first order phase transition in the 
Edwards measure of the KA model was pointed out to 
me by S. Franz. For a study of dynamical heterogeneities 
in the same model, see: S. Franz, R. Mulet and G. Parisi, 
Phys. Rev. E 65 021506 (2002). 

[18] Notice that in the inhomogeneous case - when the en- 
tropy is a functional of the density profile, the functional 
version of Eq. (7) holds, and Tsdw will only depend on 
the density profile provided its thermodynamic interpre- 
tation holds true. 



4 



